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Abstract 

An equation for the reduced density matrix which describes a free particle, that is in- 
teracting with a linearly dissipative medium, is derived using the total Hamiltonian, and 
without resorting to any artificial model. A Master equation is also obtained for the dissipa- 
tive medium. The theory erases the notion of the reservoir. It is shown that the dynamical 
interaction with the medium is not completely determined by the friction constant (or the 
corresponding fluctuation-dissipation relation). Unexpectedly, through the density-density 
correlation, the longitudinal dielectric function is shown to play a critical role. 
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The process of dissipation is almost the antithesis of unitary evolution of the wave- 
function, and therefore a quantum description of even the simplest quantum system (that 
of a free particle) interacting with a dissipative medium has turned out to be a formidable 
problem. In fact, and as is clear from recent publications - ||, there is still no agree- 
ment on the equation of motion for the reduced density matrix that describes the particle. 
The main problem has been the derivation of a suitable Hamiltonian that describes the 
whole system. So far, no one has found operators that can simultaneously a) describe 
the essential features of the dissipative medium, b) satisfy a commutation relationship, 
and c) be used to construct that part of the Hamiltonian that describes the interaction 
with the free particle. Consequently, the dissipative medium has often been modeled as 
a collection of harmonic oscillators (or similar entity ) with features prescribed by the 
fluctuation-dissipation theorem B 
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Among the many conditions on the equation for the reduced density matrix, two of the 
most difficult to fulfill have been a) an equation that is valid when quantum fluctuation is 
larger than thermal fluctuation, and b) an equation that ensures positivity of the density 
matrix. It was demonstrated that in order to describe quantum Brownian motion the 
probability diffusion coefficient typically features a time dependence, and does not vanish 
at zero temperature ||, fll2f. However those works were criticized because positivity of 



the density matrix was not guaranteed. On the other hand, in works where a term was 
added to the standard Master equation in order to ensure a positive semidefinite density 
|, @, none of the equations can really describe quantum Brownian motion 
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because the diffusion coefficient vanishes at zero temperature, and symptomatically it 
is not time- dependent. As will become clear, a key fragment of the riddle of Brownian 
motion, as well as the shortcomings of all previous approaches, originate in the certainty of 
the fluctuation-dissipation relation for the friction coefficient. In the following discourse, 
I present for the first time a description of Brownian motion that does not rely on a 
collection of harmonic oscillators (or similar) model for the lossy medium. A complete 
Hamiltonian is derived and used. 



The particular system of interest — a 
current driven RC circuit — , which is 
identical to a free particle driven by an 
external force, is shown in figure 1. In 
the absence of a resistor, the system is 
well described by the charge operator 

which satisfy 



Q, and operator $ [14 




a commutation relationship [$, Q] = 
ih. And the contribution to the total 
Hamiltonian is given by H c = Q 2 / (2(7) , 



Figure 1: Schematic of the 
circuit under consideration 
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plus the energy due to the driving force W s = —I a $. 

In this paper, it is asserted that a quantum mechanically complete description 
of the interaction with the dissipative medium can be obtained in terms of the in- 
trinsic electrical current operator, 3 d (q) = J2\rj P^iQ^x^v > with P Xr, (q) = 

2^ / d 3 re~ l <l- r {(J)* x (r)\7 4> v (r) — (f) v (r)\7 (f)* x (r)} ; and the charge density operator 
(Pcd)> Pcd( r ) = T,i eiipd( r )^d{r) = Ei ei<5(r - n). Here ^ d is the wave function of 
the dissipative medium, (f)^ the basis states, and t\ the charge of the I th specie. It can 
be shown [17]] that in wavevector space the operators satisfy a commutation relationship 
[ PcdiQ) > J d(-q)) = hqEief/m , were p cd {q) =$<PrpJr)e- i r . <n P 2> 

Construction of the Hamiltonian for the complete system started from the Hamiltonian 
for the current driven capacitor, and the interaction energies were derived using standard 
electromagnetic theory. The Hamiltonian was found to be 

H = H c + H d + W s + H je + Had + Wd . (1) 

Here H cd describes the interaction between the scalar potential Q/C and the charge 
density operator : H c d = — J d 3 rp cd (r)Q(r)/C. And the energy Hj e is due to the 
interaction between the current and field Hj e = J d 3 r Jd(r). V$(r) . The interaction 
energy of the current source driving the resistor is given by the charge density interacting 
with the equivalent potential Wd = —I s R s p c d(.q)- The Hamiltonian for the dissipative 
medium is written as: H d = £>k C k a C ka + £<? v i(Q)Pi(Q)Pcd(Q) + H <* 
, where (p^q) = J2j e 1 ^ ^) is the impurity density, Vi(q) is the potential associated with 

the scattering, and H ee is the electron-electron interaction. <np3> 

The total density matrix is written as p 2 (t) = S(t) ® p d (t) + p c (t), and it evolves 
according to the equation: ifigf- = [H, Here S(t) is the reduced density matrix 
that describes the free particle S(t) = Tr d /0 E (i), p d describes the dissipative medium 
p d = Tr s /0 E (t), and p c is the correlation between the sub-systems. In general, any 
density matrix can be written in the above form. The calculation is done in the interaction 
representation. Accordingly I write H = H + W + , with H = H c + H d , and 
W = W S + W d , and of course for any operator O, 0(t) = e% Hot Oe~% Hot . It then follows 
from Liouville's equation that 

= i W s + V.{t), S(t) ] + Tr d [W(t) + Hi{t), p c (t) ] (2) 
i% d _h = [W d + V d (t),p d (t) }+Tr s [W(t) + H, i (t),p c (t) } (3) 

i%d ~m = -[^(t)+W s ,S(t)]^p d -S® [V d (t)+Wd,p d } + [W + H i ,S®p d } (4) 
+[W + Hi, p c ] - Tr d {[W + Hi, p c ]} ®p d -S® Tr s [W + H u p c ] 
and V s (t) = Tr d {(W d (t) + H(t))p d (t)} , with V d (t) = Tr s {(W s (t) + H(t))S(t)} . 
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A detailed description of the iterative process by which these coupled equations ( |]), 
( ||) and ( H) were solved will be described elsewhere |L7|]. First a solution for p c was 
obtained (also via an iterative process), then that solution was substituted into ( 0); p d 
was computed using the improved p c , and S(t). In computing S(t), the initial values 
were obtained by assuming that initially SI and S2 are open (see fig.l), and at t=0 the 
switches are closed, then the current turned on to a value of I s . It was found that 

g = _^„ + H ,, S(t)] + U c (f + f)u c ,5, 
with U c = exp(—iH c t/h) and, 

U c § a U c h 2 = -[$, $5] Wj + [$, QS] - [Q, QS] W+-^ 

+ [$, Sd>] Wj - [$, 5Q] Wjii + [Q, SQ] W-± 

U c S b U c h 2 = -([$, [Q, S}} + [Q, [$, ^])~Re(W£) - ([Q, {$, 5}] - [$, {Q, S}])~Im(W£) 

The spectral densities of operators of the dissipative medium, and functions 
related to the spectral densities appear in equation ( [|). They are defined 
as follows: W}(T,t) = ^Tr d {p d (r)4(q,r)J d (q,0)}L 2 dr ; W+(T,t) = 
^rTr d {p d (T)4(q,T)J d (q,0)}L 2 dr ; W;(T,t) = ^Tr d {p d (r)p cd (q,T)p cd (-q,0)} dr 

Other integrals of ensemble averages are similarly defined: Wj(T, t) = Jq < 
4(q,0)J d (q,r) > dr Lf ; Wj t (T, t) = / *r < 0) J d (q, r) > dr ; and 

W~(T, t) = Jq < p cd (—q,0)p cd (q,r) > dr. The integral of the cross-correlation function 
is given by W+j(T, t) = J* < p cd (q, r)J d {q, 0) > dr. 

In the second step, the solution for S(t) is used in the computation of p d {t). Hence 
now V d (t) = Q d Tr s (QS(t))/C + {J d - I s )Tr s {<f>S{t)), and from ( |) it follows that 

ih ^ = [v d (t)+W d (t),p d (t)] + ^J*dt?T^^ ] 

- i j[V Tr.[#(f),5(0 ® fa(0, AiC*)] ] • (6) 

The spectral density of J d and related functions were evaluated from the relationship 
between the conductivity and the retarded current-current correlation ( II a ^) function 
0: 

Re a a/3 (q,uj) = - Im U a p(q, lu)/uj (7) 
n a/3 ( 9 , w ) = -_ / dte{t-t')e^-^ <^ d \ [4 a (q,t),J d/3 (q,t')} \^ d > 

V i J-oo 

where V Q \ = Ljj is the volume. The trace was taken with (p d (t) = p d (0)), and the 
ensemble averages were taken in the long wavelength limit (q — > 0). 

The spectral density of the charge fluctuation {W^ 1 ) in the resistor plays a critical role 
in determining the off diagonal elements of the density matrix. This charge fluctuation 



4 



is the fluctuation of the conjugate of the intrinsic current operator, and is not the charge 
fluctuation on the capacitor. The latter is determined mainly by the current fluctua- 
tion. Now the density-density correlation function is related to the longitudinal dielectric 
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(hS + i(QS$ - ^SQ))^ \m{W+) . (9) 



, , 1 - . lim — — / dre^ T < T T pjq,r) pj-q,0) > . (8) 
e\\(q,ou) iu» -» u+i5 V i \q\ Jo 

Evaluation of W*(T,t) would have been straightforward if the longitudinal f-sum rule ( 
linig^o Im e~ l (q, oS) = — ttu p [5(u> — u> p ) — S(u + u> p )}/2) could have been utilized. But 
such is not the case since the low frequency variation of the integrand is important (see 
also ( [TO]) below ). In the high temperature or semiclassical (h/3uj << 1 ) limit we find 
lim^o < Pcd(q,T) p cd (-q,0) > = ^e"^™ /(4R s r m ) , and hence W+(T,t) = h 2 (3(l- 

e~*/ Tm )/ (4i2 s ). Here r m is the momentum relaxation time of the electron. <np4> 

Hence, the Master equation for S(t) can be written as 

f)Q 1 1 1 j 

a = a [ 2C Q2 ~ Is ^ s{t)] ~ [$ ' [ *> s]] ¥ Re {wj) + ^'W'^^M^) 
" [g ' [g ' s]] Jfky Re{w ^ + [^[Q,s]}-^Re(w+) 

- [$,{$,5}] ^ lm(W]-) - [Q,{Q,S}\ j^y 2 Im(H/;) 

2 

iFc 

It is worth noting that this equation is very clean: the coefficients of the various terms 
do not depend on (nor were they chosen to match) any classical/quantum criteria that is 
outside the scope of the Hamiltonian, and furthermore the coefficients are directly related 
to Green's functions of the current and charge density operators. There are more terms 
than expected. But more importantly, the physics contained in those terms is clear and 
fundamental. The terms proportional to Im(W^), and Im (Wj") yield the change in 
the density matrix due to the change in the rest mass of the particle. Examination of 
Im (Wj) shows that it is similar to the calculation of the renormalized mass of a free 
electron that is interacting with a radiation field. 

The term L a = [Q,[Q,S])qzl , has received a lot of attention, since supposedly it 
is required to preserve positivity of the density matrix [f[5fl , |I6| , [|J, 0. And despite 
significant efforts (JH, {§, [@]), the coefficient of that term has only been obtained at high 
temperature (at, = oiLi/iksT), thereby leaving the question of positivity still open at 
low temperatures. Furthermore the constant an is not the same in the different works. 
It was shown by Lindblad that positivity of p(t) is assured for Master equations of the 
Lax-Louisell form: p = —ih^lH, p\ + X^([Aj-, pAj\ + H.C) (where Aj is composed 
of system operators). Under the conditions implicitly considered by the aforementioned 
investigators (high temperature and the dissipative medium is not quantized), equation 
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( is of the Lax-Louisell form. Under all conditions equation ( ||) is of the Lax-Louisell 
form, and ( |9|) would have been of that form if terms proportional to the real part of 
the cross correlation were kept. However when the correlation functions are computed 
for the dissipative medium in equilibrium, then ReWlj = 0. If the dissipative medium 
is driven very far from equilibrium then the aforementioned terms can be reinserted, 
and the ReWw computed using equation ( |1^). The really astonishing result concerning 
L a that is reported here is that the coefficient {oll) is related to an important physical 
quantity: the longitudinal dielectric function (or the density susceptibility function when 
the system operators are position and momentum). Hence with ( |8|) it is now possible to 
give an exact expression for (a^), that is valid at all temperatures: 

M Kit) ) = - » v„ Jim r Im ( _L_) dw . (10) 



(2tt) 2 g-o J-oo 1 - e-f 3 ^ \e\\(q,u)J lu 

The physical mechanism associated with L a is the screening of coherent charge fluctu- 
ation (on the capacitor) by the medium. One can say that the friction (R(u)) damps 
the ensemble- averaged charge/momentum, while the coherent momentum fluctuation is 
damped by the density susceptibility. 

In a semiclassical computation, the diagonal elements <u\S\u> can be treated as a 
quasiprobability distribution, here \u > is the eigenvector of the charge Q\u >= u\u >. 
The term L a does not affect evolution of the diagonal elements. From ( P), the evolution 
of < ulS'l'U > is given by a Fokker-Planck equation, with a probability diffusion coefficient 
D uu = Re (Wj). In the purely ohmic limit (r m — > 5), D uu ~ K B T/ R s + h/ (7rR s t), and 
hence for all times t < h/3/n, the diffusion and particle's evolution are purely quantum 
mechanical. Evolution of the ensemble average < Q >= Tr{QS(t)}, is given by 

i±§> - ^ lm(wll) <q>=i, (n, 

which reduces to the classical equation: Q + (R S C)~ 1 Q = I s . Equation ( [11]) represents a 
true milestone: the derivation — for the first time — of the classical equation that describe 
a dissipative system from a purely quantum theory. Note that Wjl is directly related 
to ipd, the wavefunction of the dissipative medium. In theories that use a collection of 
harmonic oscillators, the classical equation is not derived, but rather it is used to fit the 

free parameters (spectrum, etc) of the theory. <np5> 

It is now possible to write — for the first time — a reduced density matrix that describes 
the evolution of a dissipative medium that is coupled to a quantum object. Continuing 
from ( |B|), I obtain: 

^ = l[H d + W d + V d ,p d } - [l d Ah,P d ]]^eWUt) 

+ ( [Q d , {I d , p d }} - [I d , {Q d , p d }} ) ^Im Wfc(t) - [Q d , [Q d ,p d }]± W+ Q {t) 
+ [l d ,{I d ,p d }\^lvxWUt) . (12) 
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The charge operator associated with the charge density is defined as: Qd = 
\imq^ p cd (q), and Id (amps) is the long wavelength limit of the current through the 
resistor. In this bistratum quantization approach, one set of operators describe the in- 
dividual modes or individual particles of the dissipative medium (C' x s, see numbered 
paragraph 2 (or< np2 >)), while the key operators that describe the coupling and enter 
the Master equations 1) satisfy a commutation relationship in wavevector space, and 2) 
are composed of the microscopic operators. Equation ( [12]) still contains the full Hd 
and hence a direct solution is far from trivial. Nonetheless, the equation is manageable 
because the results of the equilibrium (only Hd) solution are known. The partial spectral 
functions are given by: 

W+(T,t)= f dt'Tr^S^K^t-t')^) - Tr s (S(t')A,)Tr s (S(t')K(t -t')) 
Jo 

with Ai = $, and A 2 = Q/C . 

Equations ( ||) and ( [12]) radiate the sublime grace of the new method: bistratum 
quantization, in concert with a coupled systems approach that recognizes both objects as 
quantum entities. The foundation of theories of dissipative quantum systems has so far 
been inscrutable. In contrast, here the quantization of the medium, quantization of the 
operators used to construct the interaction Hamiltonian, and the fundamental connection 
between those operators and the microscopic operators are all clear. Parenthetically, it 
is possible to draw a distinction between an open system (typically found in quantum 
optics), and an inherently dissipative system. In an open system such as an atom in a cav- 
ity, a process such as spontaneous emission is sometimes viewed as dissipative, but if the 
number of modes is reduced then the process becomes reversible. In contrast, resistance 
to even dc current flow, and irreversibility is present in the typical closed system. The 
difference is due to the underlying mechanisms in the dissipative medium: in one case a 
set of independent modes, and in the other a strongly coupled (H ee ) electron system. It 
is important to note that an intrinsically dissipative medium is not really quantized in 
models that employ a collection of oscillators (or similar). One may think that under the 
quantization of the harmonic oscillators, ... but quantization in a space is irreducible, 
there is nothing under, .... under. A direct consequence being that the evolution of the 
dissipative medium could not be determined, and hence questions concerning operators 
of that medium were part of the unspeakable. 

In summary, equations for the reduced density matrices of two interacting sub-systems, 
a free particle and a linearly dissipative medium, were derived. The approach to the 
40 year conundrum (see <npl> ) of a Hamiltonian description of the interacting sub- 
systems, and quantizing the whole system is clear. Operators that satisfy a commutation 
relation — in wavevector space — were utilized to construct an interaction Hamiltonian. 
The total density matrix was written as a tensor product of the reduced density matrices 
plus a correlation matrix. Partial spectral densities that appear in the Master equation 
for the particle are given in terms of well known quantities: the retarded current-current 
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correlation function, and the retarded density- density correlation function. 

The theory is particularly important because a) in many cases of interest, operators 
belonging to the dissipative medium carry the required information and now questions 
that contain those operators can be answered, and b) it appears to be the basis of a 
truly quantum circuit theory. For example, in the system shown in figure 1 (which may 
be part of a larger circuit), the mean current through the resistor, and the fluctuation 
of that current often correspond to the detected signal and the fluctuation of the signal. 
A circuit theory that can describe quantum transport is particularly important these 
days, because everybody's vision of the future now includes nanoscale and macroscopic 
quantum circuits. 

One of the most surprising features of Brownian motion that was unveiled in this 
investigation is the role played by the longitudinal dielectric function. It is easy to un- 
derstand why all previous investigations have missed this important feature. Classical 
and semiclassical studies are based primarily on the Langevin, or quantum Langevin 
equation (for the momentum), and as can be ascertained from <np5>, the density cor- 
relation does not contribute to the Langevin noise force. Quantum studies were based 
on reproducing the friction constant R(u), either directly through the Caldeira-Leggett 
model, or indirectly through a fluctuating force and the fluctuation-dissipation theorem. 
However as we have seen the dynamics of the medium is underdetermined by R(u>), or the 
associated fluctuation-dissipation theorem. How reliable are seemingly exact quantum 
models of an object, where exactness is determined by some classical criteria? It appears 
that a quantum object can only be known via the correct Hamiltonian. 
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